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ABSTRACT 

We present QCD corrections to the Drell-Yan process in the transversely polarized, longi- 
tudinally polarized and unpolarized cases. The analytical results are presented in a form 
valid for all n-dimensional regularization schemes. A universal mass factorization scheme 
is presented in which the results reduce to those of dimensional reduction. The connection 
between the parton distributions and fragmentation functions of dimensional reduction and 
those of dimensional regularization is elucidated in a simple manner. Numerical results are 
presented for proton-proton collisions at energies relevant to RHIC (Relativistic Heavy Ion 
Collider). The perturbative stability of the transverse and longitudinal asymmetries is 
investigated. 

PACS numbers: 12.38.Bx, 13.75. Cs, 13.85. Qk, 13.88. +e 
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I. INTRODUCTION 



The unpolarized Drell-Yan process has been studied rather extensively in the hterature, 
including 0{as) |l|, and 0{a^) corrections. As well, the 0{as) corrections to the 
corresponding longitudinally polarized |Q and transversely polarized processes ||^, Q have 
been studied. What was still lacking is a unified picture for dealing with the polarized 
processes. The basic problem is the ambiguity associated with defining the 75 matrix, 
or s^'^^'' tensor, in n dimensions; both of these objects arise in polarized processes. For 
unpolarized QCD processes, dimensional regularization (DREG) preserves all the necessary 
invariances and symmetries to do calculations to any order in Hence DREG is the most 
commonly used regularization for QCD. The ambiguity associated with the continuation 
of the 75 matrix makes it impossible to uniquely define higher order corrections (HOC) 
for polarized processes using DREG. Various prescriptions are available, but problems 
with either mathematical or physical consistency generally arise. As a result, another n- 
dimensional scheme, dimensional reduction (DRED) may be used. This scheme avoids the 
75 problem, although it requires certain ultraviolet (UV) counterterms which are the same 
in both unpolarized and polarized processes and may be unambiguously determined. 

In this paper, we present analytical results for unpolarized and (both longitudinally and 
transversely) polarized Drell-Yan in a form valid for all ra-dimensional schemes. For the 
polarized case, the ambiguity (or scheme dependence) in the DREG results is parameter- 
ized by the ambiguity in the polarized n-dimensional split functions. As well, we present 
numerical results for p-p collisions relevant to the Relativistic Heavy Ion Collider (RHIC). 

We go further to show that for a wide class of subprocesses, including the one-loop 
corrections to Drell-Yan and deep-inelastic scattering, DRED is simply equivalent to a 
particular mass factorization scheme in DREG. We call this scheme MS^ (or MS^) since 
it involves subtracting the e-dimensional part of the n-dimensional Altarelli-Parisi split 
functions, where n = A — 2e. As a consequence, the final results in the MS^ scheme are 
regularization scheme independent within the n-dimensional schemes. The final result is 
equivalent to that obtained in DRED and all ambiguities associated with the continuation 
of the 75 matrix are subtracted via the n-dimensional split functions. 

We will also show the connection between the DRED parton distributions and frag- 
mentation functions and those of DREG in a simple manner. More specifically, we show 
how to convert existing DREG distributions into ones suitable for use with cross sections 
determined using DRED. This is important since DRED is equivalent to 4-dimensional he- 
licity amplitude techniques which considerably simplify perturbative calculations. We may 
thus calculate new unpolarized cross sections using DRED or helicity amplitudes and then 
simply convolute them with the DRED distributions obtained from well-known unpolarized 
DREG parton distributions and fragment at aion functions. 

Similar conclusions (for unpolarized processes) may be obtained in the approach of |]^, 
which converts DRED cross sections into DREG ones by considering differences in the 
Lagrangians and using fictitious £-scalars to calculate the differences in the cross sections. 



1 



Transition rules between the two schemes are also given in [^]. Here, we take a simpler 
and more phenomenological approach, investigating how the scheme dependences arise 
in the Feynman graphs. The connection between the two schemes is simply the relation 
between the distributions of the respective schemes. This allows for easy interpretation and 
extension to a wide class of processes. We also explicitly consider polarized observables, 
unhke 0, 1. 



II. n-DIMENSIONAL REGULARIZATION SCHEMES 

There are two parts to the dimensional continuation: the continuation of the momenta 
and the continuation of all other tensor structures (i.e. gamma matrices). The continuation 
of the momenta is unique, but there are various methods for continuing the tensors. The 
choice of the latter defines which dimensional method is being used. 
Continuation of the Momenta 

For the continuation of the momenta, all momenta and phase spaces are continued 
to n dimensions ||^, |10[- The phase space integrals are continued by generalizing integer 



dimensional integrals to non-integer dimensions. Consequently, all loop integrals can be 
reduced, using Feynman parameters, to the fundamental integral 

"^"^ ^^'^^ - ■C"-™+"/2^(r + n/2,m-r-n/2), (1) 



(27r)" (g2 _ C)™ (47r)'^/2r(n/2) 

(see, for example [|ll|]) with m > 0, r > and B the Euler beta function. As well, defining 
n = 4 — 26 {n' = A ~ 2e') with e < {e' > 0) we see that e' is required for UV divergent 
integrals and e for infrared (IR) divergent ones, initially. Then we must continue to 

(2) 

since we can only work in one dimension at one time. From (||) it follows that massless self- 
energy insertions on massless external lines vanish. This means that on-shell wavefunction 
renormalization is trivial when all particles are massless. Hence, in the absence of coupling 
renormalization (i.e. gluon self-energies), effectively no UV renormalization is required. 
Continuation of the Tensors (DREC) 

In DREG, one continues the metric tensor and the gamma matrices to n dimensions. 
Letting gl^^ denote the ra-dimensional metric tensor, we have the relations 

g'ng";. = + = 2^7r- (3) 

As well, the usual convention is to take 2 — 2£: = n — 2 helicity states when averaging over 
initial gluons/photons. This is related to the continuation of the helicity sum rule 

E^^(p,A)A-(p,A)--^?r. (4) 
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Here A^{p, A) is the gluon/photon polarization vector for gluon/photon momentum p and 
helicity A. Different conventions simply amount to finite renormalizations of the parton 
distributions (which we will see arise from differences in the n-dimensional split functions). 

There exist two popular methods for continuing the 75 matrix (e'^'^^^) tensor to n dimen- 
sions: the anticommuting-75 scheme |]12[ (see also |13| concerning e'^'^^'') and the HVBM 
scheme [|], . 



In the anticommuting-75 scheme, we use the relations 

757m = -7/.75, 75 = -1- (5) 

If traces with only one 75 occur though, there are known mathematical inconsistencies [0]. 

In the HVBM scheme, we formally take n > 4 (with regards to the tensor algebra) and 
keep the 75 and e'^'^'^p in 4 dimensions so that 

{75, 7/.} = 0: ;U < 4, [75, 7^] = 0: /i > 4, (6) 

which follows from the definition 

75 = (7) 

where 

^Ml/12Ai3M4 f^i ^ 4, (8) 

otherwise, it is the usual Levi-Civita tensor. 

This scheme is mathematically consistent, but cumbersome. Physically, it has the 
problem that the non-anticommuting 75 leads to non-conservation of helicity of massless 
fermions in a minimal subtraction scheme like MS |TH . 
Dimensional Reduction 



Dimensional reduction \W\ is perhaps the simplest of all the dimensional methods. It 
was originally introduced because DREG violates supersymmetry. As will become obvious, 
it is also manifestly mathematically consistent. The idea is simple; all 7-matrices and ten- 
sors are taken to be 4-dimensional, and formally n < 4. This implies that the components 
of all momenta between n and 4 must vanish. We have the following contraction identities 

g^'^g,, = 4, gil'^g;, = g^'^g;, = n (9) 

and the usual 4-dimensional relations like 

7V + 7V = 2^''"- (10) 

It is also useful to define 

^^^^^^[g^^ -g^^). (11) 

This method is particularly simple for the calculation of tree graphs (i.e. graphs not 
involving loops) since the traces are equal to their 4-dimensional counterparts, implying 
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gauge invariance. One may thus use 4-dimensional helicity amplitude methods, for instance. 
Then the phase space integrals are carried out in n dimensions, providing and IR regulator. 
As well, the anticommuting 75 implies helicity conservation of massless fermions. 

The only subtlety comes from the fact that the virtual momentum integrations generate 
the tensor g!^'^ , which is generally contracted with 4-dimensional 7- matrices. This can lead 
to a term ~ 7^ which gives the incorrect Lorentz structure and must be removed by 
a counterterm. In fl^ the counterterm for the quark-7(Z) vertex was presented. It is 



(working in the Feynman gauge) 



with Cp = 4/3 (i.e. the Feynman rule for the counterterm is obtained by making the above 
substitution in the usual rule). For the lepton-7(Z) vertex, we use (^2|) with 

Cpg^ - e\ (13) 

Throughout, we consider (p!2[) to be a Feynman rule for DRED. For the type of processes 



considered here, ([l2|) is the only counterterm required to make DRED physically consistent. 



III. ANALYTICAL RESULTS FOR DRELL-YAN 



We will first consider the unpolarized and longitudinally polarized cases, then the trans- 
versely polarized case. We have the general process 

A(Pi, \a) + 5(P2, Ab) ^ l-{pz) + /+(P4) + X, (14) 

where A^, A^ denote the helicities of hadrons A, B. The unpolarized and longitudinally 
polarized cross sections are defined, respectively, by 

a = i[a(+, +) + a(+, -)], Aa = +) - a(+, -)] (15) 

in the notation cr(A^, A^). 



The general 2^2 [2^3] subprocess contributing to (14) has the form 



a(pi, Ai) + 6(P2, A2) ^ 7*(g) + [cik)] ^ /-(pa) + l^ip^) + [cik)] (16) 

for general partons a, b, c. 

Firstly, we define the process-level invariants 

S=iPi + P2)\ = (p3 + pa)', r = M'/S. (17) 
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In the parton model, we have 

Pl = XaPl, P2 = XbP2- (18) 

Hence we may define the subprocess invariants 

^2 n W T 
S = {P1+P2) = XaXbS, W = = = . (19) 

S oXa^b ^a^b 

The unpolarized [polarized] momentum distributions are given by 

[A]F,/i{x„ M}) = x,[A]U/i{xi, M]) (20) 

where the [Aj/j// are the unpolarized [polarized] parton densities for parton i in hadron /, 
evaluated at factorization energy scale M|. 

The parton model expression for the Drell-Yan cross section corresponding to (n) is 
then 

IM^ - 5/;^£ t^|Alf./.(^.Mj)[AlF,.(..M|)|Al^ (21) 

where 

Wi = T/Xa, Xb = Wi/w (22) 

and [A]o"af, is the unpolarized [polarized] subprocess cross section corresponding to (|l^). 
We must consider the subprocesses a = q,b = q, c = g and a = q, b = g, c = q, which are 
symmetric under a ^ b and g ^ g as far as [Aldaab/dM"^ is concerned. 
Define the unintegrated leptonic tensor as 

<RED = f^''e'[p^p^,+p2p§-{My2)g-^] 

<REG = f^'^e'[p^p^,+p2p',-{My2)gf] (23) 

where the arbitrary mass scale /i^^ arises from the n-dimensional coupling e^/i^^. 
Furthermore, define the integrated leptonic tensor as 

~ 7 (27r)-i2p3,o M2 ^ ■ 

One finds 

- ,2. 2^ Mg^)'^ r(i-e) 1 

^DRED e 2^ ^2 r(l-2£)(3-2£)(l-2£) 

(25) 

The corresponding DREG tensor is obtained by replacing (7"^ — > g^^. This gives 



[l-ef-^ + \-{^-2ef— 



..p _ ,2. 2^^(?')"' r(l-e) (1-e) 



''-DREG — p e 



2" 7r2 r(l-2e)(3-2e)(l-2e) 



a 3 

q2 yn 



(26) 
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The part ~ does not contribute to the cross section, as can be seen from gauge 

invariance. Hence the integrated leptonic tensor is effectively a constant. Nonetheless, we 
keep all the terms for completeness. 

We may then define the unpolarized [polarized] subprocess hadronic tensor [AJll/"^^ 
through the unpolarized [polarized] subprocess squared Feynman amplitude 

[A]|M|^, = ^L.,[A]W^-^ (27) 

where A|M|^j are defined analogously to (plSf). Having done so, we may write the 

2 — >^ 2 phase space as 



Similarly, for the 2 — * 3 phase space. 



r , -| ^g"afc,2^3 _ 1 



TT 21-2^ T(l-e) 



j dyy-\l - yr[^\W:,^^,^,C^p 
^ 



(29) 



where ?/ = (1 + cos^)/2 and Q is the angle between pi and k in the pi, p2 cm.. This is all 
we need to calculate [Ajficr^^/dM^. 

In order to present [Ajrfaab/dM^ in a form valid for all n-dimensional schemes, we 
must first give the general form of the n-dimensional split functions Pjj(z), related to the 
probability of parton j splitting into a collinear parton i having momentum fraction plus 
an arbitrary final state carrying the rest of the momentum. 

We may write 

PlAz, e) = P,^{z, e) + 5(1 - z)PUe), (30) 



with 
and 



P.^i^,e)=P<'\z)+eP<'%z) (31) 
P^is) = p!-' + sP^\ (32) 



In DRED, P,f'{z) and P-f are zero. In other words 

i^°^^^(z) = /^5(.), (33) 

where P^j{z) is the usual 4-dimensional split function. 

One might wonder how to determine the P^{z, e). It is done in the same way as for the 
4-dimensional case, but keeping the terms of 0{e). 

We can make this clearer by considering the 2 m process (all particles massless) 

a(Pa) + KPb) Ci{ki) + C2{k2) H h Cm{km)- (34) 
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When ki is coUinear with pa, we have (at one- loop) 

[A]\M\l_,^{ab^c,---Cm: 
where a — > d + ci and 



d Pa - 1^1 



fcl ^ (1 - Z)pa ^Pd^ ZPa, 



Z < 1, 



(35) 



(36) 



with APjj being the corresponding polarized split function. Then P^j is determined us- 
ing probability and momentum conservation, and it only appears when i = j. Also, by 



definition APf- 

''J 



In this paper, we will need [A\Pqq{z) and [/^Pqg{z). In 4 dimensions flS 

2 



1 + '^" 3^,^ 



1-2 



1 



'1 - 2z + 2^' 



2V /' <79V / 2 

For the unpolarized case in DREG, the e-dimensional parts are unique [11^ 

2 



(37) 



-Cf(1 



'J? 



z"^ — z. 



p: 



19 



0. 



(3^ 



except that P^g'^ depends on the convention used for averaging over initial gluon states. In 
the anticommuting 75 scheme, one has 



= AC.(^) = K.,S^) - PliM 



<M 



1+1+^ 



(39) 



(with the ±'s denoting helicities) as a consequence of helicity conservation (i.e. an anti- 
commuting 75). This is not true in the HVBM scheme due to (|^), which violates helicity 
conservation of massless fermions (see [|l^] concerning the polarized split functions in the 
HVBM scheme). We will show that this problem is overcome in the MS^ factorization 
scheme, which we shall now define. 

Factorization of the mass singularities is equivalent to expressing the bare parton distri- 
butions (and fragmentation functions) in terms of the renormalized ones. In the MS (MS^) 
scheme, this is done via 



where 



c{e) 



li/A 

1 /47r/i 
1 {Mf 



E ^ r-m7/A \y^ Mj)[A]Pf (a:/y), (40) 



2\ e 



r(i-e) 
r(i-25) 



^1 (--7E + ln47r) + 0(.) 



(41) 



The conventional factor [ji^ /M'jY is not necessary, but it allows for a distiction between 
coupling renormalization and mass factorization energy scales. We will take Mf = fi in 
our calculations. In other processes though, this distinction might be necessary in order to 
avoid large logarithms. 

For the fragmentation functions Va/i, representing the probability for quark i to split 
into hadron A, the corresponding renormalization is 

[AMdz) = [A]vYjt\z,M}) + ['^[A]V%^^\y,M])[A]Pj:^{z/y), (42) 

c j ZTT Jz y 

where aj = a^, unless j {=A) = 7, in which case aj = a. 

It is clear that the MSg scheme is just the MS scheme, extended so as to subtract off 
the entire n-dimensional split function. In DRED, MS^ is equivalent to MS since there is 
no e-dimensional part of the split function. 

We are now in a position to write down the results for [A]daab/dM'^ in a form valid for 
all n-dimensional schemes. We start with the qq subprocess. The Born term is given by 

[A]^^ ^ [A]xB{e)S{l - w) (43) 

with 

DRED/ X _ _AvDRED/ N _ 6^ 27r^+" 4^ (2 - g) r(l - g) 

Xb [e)- ixxb KS) ^^2-2-^4+2./^ (3-2g)(l-2g)r(l-2g) ^ ^ 

and 

xT^'^ie) = 2^^xi^"°(^)- (45) 
Here Nc — 3 and Cq is the quark fractional charge. In the anticommuting 75 scheme 

^X^%) = -xT^'^ie). (46) 

This is not true in the HVBM scheme though. Of course, in the limit e ^ helicity 
conservation is restored in all schemes so long as there are no 1/e poles arising from mass 
singularities. If there are such 1/e poles, then one needs a scheme such as MS^, as we will 
see. 

The factorization counterterm in the MS (MSg) scheme is 

lAl^ ^ H|A]x.(.)S»-nAie'(») ( ^ (4.)'- (47) 



'ciM2 ^^""^ ^27r ^ j <jg v ; ^ ^2 y r(l - 2£) {Mj 

The gluonic bremsstrahlung contribution is 



X 



l^Sd - .) - HI^SM + 8 f i^i^) - 4(1 + - - 2(1 - 



The virtual contribution is 



Til 



2n \ I r(l - 2e) 



2 



(49) 



So, adding ( ^3l) and (|47|) - ( ^9|) we obtain the total result 



as 



[A]xBiO)[Sil-w) + CF—w 



7\5{l-w) 



(50) 



+ 



/ln(l 



w 



\ 1 



w 



(w) S 



w 



2(1 



w] 



[A]d. 



qq 



where 



MSe 

qq 



0, 



Op 



(51) 



We see that the 0{e) scheme dependence of the Born term cancels with the 1/e (1/e^) 
divergences multiplying it. Also, we see that all the n-dimensional regularization schemes 
give the same answer in the MS^ scheme and it corresponds to the DRED MS answer. 

We now consider [A]dcrqg/dM^. There is no C(l) term (in a^). At 0{as) there is a 
factorization counterterm contribution which is given in the MS (MS^) scheme, by 



d(f^^ 



\[A]xB{e)^w'^^'[A]P^-\w) 



2\ E 



rfi 



r(i -2£) \M] 



(52) 



The bremsstrahlung contribution is 



a. 



'AlTfl 



2\ ^ 



ra-e] 



M2 / T{l-2e) 



X 



--[A]P<{w) + 2[A]P^^iw) ln(l - w) 



w 



+ w(l — w) 



(53) 



Adding 

da. 



and 



we obtain the total result 



[A]- 



qg 



'dM^ 
where 



[A]XB(0)^t. [A]P;(^) 



liiT72+21n(l-^) 



Mi 



w 



[A]d, 



qg 



0, 



MS 

qg 



-[A]P. 



qg ' 



w 



:i + 3w) + [A]d,,^, 
(54) 

(55) 



In both the qq and qg cases, we verify that the unpolarized DREG MS result agrees 
exactly with that previously determined (see for example 0). Since 7^ AP^ in the 
HVBM scheme, we see from (|50|), (|51|) that ddqq/dM^ ^ -Adagg/dM^ in the HVBM 
scheme if one uses MS. But this is a physical requirement. Hence, if one uses HVBM 
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regularization, then it is necessary to use a subtraction scheme hke MS^ or one which 
subtracts at least the hehcity non-conserving part, AP^''^ — P^'^ as well (see also [0) in 
the polarized case. Of course, it makes more sense to subtract the entire AP^^ since this 
leads to regularization scheme independent results. 

Now let us consider the Drell-Yan process with transversely polarized hadrons {trans- 
verse Drell-Yan). The general process is 

A(Pi, Si) + 5(P2, ±^2) ^ l-{Pz) + 1^{P^) + X (56) 

where the Si are reference spin vectors satisfying 

= -1, = ^,J = 1,2 (57) 

implying that 5*1 and 5*2 lie in the plane transverse to the beam axis. Now, letting | denote 
polarization in the direction of the spin axis and letting | denote polarization opposite to 
the spin axis, we may define the transversely polarized cross section as 

ATa = i[(T(T,T)-a(T,i)], (58) 

in the notation a{Si, ±S'2). 

The general 2 — 2 [2^3] subprocess contributing to ( p^ is 

g(Pi, si) + g(p2, ±S2) ^ l*{q) + W)] ^ /-(ps) + /+(P4) + W)] (59) 

with 

si = Si, S2 = S2 (60) 



and At(3" defined analogously to (pSl). There is no qg subprocess since gluons cannot be 
transversely polarized. We may define the transversity distribution 

ATfg/Aix,M]) = AtF,m(x,MJ)/x = f,^/A,ix,M]) - f,jA,ix,Mf), (61) 

which is often denoted as hl{x,Mj). 

Let 03 denote the azimuthal angle of about the beam axis (with respect to some 
reference axis) and let 63 denote the angle between pi and p^ in the pi, p2 center-of- 
momentum frame. Then, the quantity we are interested to calculate is A^da / dM'^dcf)^,- 
This was done in using DRED, where A^^da / dM'^dVt^ was first determined, with ^3 
representing the solid angle of p^ in the cm. of pi, p2- Then A^^da / dM"^ dcf)^ was obtained 
via 

= / (2 COSt^3 —. 62 
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Of course, if one was only interested in A^da / dM dcj)^,, then one could integrate over 6, 
separately for the Born term, loops, bremsstrahlung and factorization counterterm, then 
add all the different parts to get a finite result for ATdcr/dM'^d^)^. Either way, the expression 
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for ^•j:(J AB / dM'^ d(f)^ is obtained from (|2T| ) by replacing [A] At and differentiating with 
respect to 03. 

From the form of the unpolarized and fongitudinally polarized results, it is straightfor- 
ward just to take the final DRED result of and put it in a form valid for all regularization 
schemes. The result is (with 0i, 02 the azimuthal angles of si, S2) 



AT^da 



a. 



27r2 



- 7 



/ln(l 



w 



\ 1 — w 



+ 2 



A^Ptiw) 



qq\ 



|„__81n(l-u.) 



6w- 



In w 
1 — w 



+ 4(1 - + Axrf 



with 



and 



AtXB = COS(01 + 02 - 203 



(ATP</ + P,t^(l-t.)) 



<?<? 



(63) 



(64) 



(65) 



The transversity split function A^P^^ is obtained via (|35|) with [A] ^ Ax. In 4 dimensions, 
AxPi, is given by 



AtPHz 



qq\ 



2 + -5{l-z) 



In the anticommuting 75 scheme, it is straightforward to obtain 

AtP</W = -Cf(1-^). 



(66) 



(67) 



The transversity renormalizations corresponding to (0), (^4 ) are obtained by replacing 
[A] — > At. We verified explicitly that the form of (|40| ) does indeed hold for the transverse 



case. 



IV. CONNECTION BETWEEN n-DIMENSIONAL SCHEMES 



Here we will show how to convert results of one n-dimensional scheme to those of 
another in a straightforward manner. We do this by examining the origin of the scheme 
dependent parts. Strictly speaking, this only applies to processes not requiring coupling 
constant renormalization, since other differences may arise from the UV sector. On the 
other hand, it has been shown ETI, |22[ that the UV sectors of DREG and DRED in QCD 



can be related via a finite 0{a1) renormalization of the coupling. Namely, 



^DRED ^ ^MS 



27r 6 ' 



(6^ 
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(see for example [Q). In other words, one may go from one scheme to the other by simply 
expressing the coupling of one scheme in terms of that in the other. We will assume this has 
been done so that the only differences may arise from the IR sector. Then all the following 
argumentation can be seen to apply to all one-loop QCD processes. 

As an example, we will consider the qq subprocess in unpolarized Drell-Yan, to show 
the origin of the scheme dependences. Then we will show that the same argumentation 
holds for all one-loop processes. 

In order to extract the scheme dependent parts, we need only consider terms which give 
rise to 1/e poles. This is because the scheme dependences come from ^ ■ 0{e) terms, where 
the 0{e) terms are in general scheme dependent. We therefore consider the contribution 
to dcTqq/dM'^ whcu k is coUinear with one of the initial partons, say pi. From (|35|) we see 
that 

12 ^5l(^'^)|,^|2 



|M|^^3~^i^|M|^„,„ w<l (69) 
Pi ■ 1^ 



with 



k^{l-w)pi. (70) 
After phase space integrations, this will yield a contribution to daqq/dM"^, 



e \ 2e 

Hence (noting that P^''^{w) ~ 1—w) we get a scheme dependent part from the bremsstrahlung, 

dff^^ 

^-xsmn^) (72) 

as well as an equal term arising from k ^ {l—w)p2- The soft divergent terms ~ (^5(1— w))-| 
cancel exactly with the loops, having the same overall factor, and hence do not lead to 
scheme dependences. 

As well, there is a scheme dependent term coming from the loops. By definition of Pqq, 
there must be a term proportional to 

This will lead to a scheme dependent part 

^ - XB{0)P'q^S{l-w). (74) 



dM^ 



Hence, the entire scheme dependence can be traced back to the process-independent 
n-dimensional split functions (their e- dimensional part). We also see explicitly why the 
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MSe factorization scheme will lead to regularization scheme independent results; all the 
regularization scheme dependent parts are subtracted. Of course, if one has longitudinal 
or transverse polarization, all the above holds with 

XBi^) A[t]Xb(£:), Pqqiw,£) A[T]Pgg(w,e). (75) 

For conciseness, we will drop the A[t]'s with the understanding that the same argumenta- 
tion holds for the polarized cases. 

It is now clear how to convert results calculated in DRED to the corresponding DREG 
MS results, if desired. One simply replaces (defining Pfj{z) = P^''^{z) + Pj^'^5(l — z)) 

Piiz) ^ Piiz) - sP!,{z) (76) 

in the factorization counterterm. 

Or, to convert from DREG to DRED (i.e. MS^), simply replace 

Piiz)^Pi;{z,e) = Piiz)+eP!^iz) (77) 

in the factorization counterterm. Hence one can go from any n- dimensional scheme to any 
another using (|7|) and From the previous argumentation, it is clear that polarization 
poses no difficulties in this approach. This procedure is equivalent to expressing the parton 
distributions and fragmentation functions of one scheme in terms of those in the other, as 
will be shown in the next section. 

Up until now, we have used Drell-Yan as an example. We now show that the argumen- 
tation here applies to all QCD processes at one-loop once the UV sectors have been made to 
agree (if coupling constant renormalization is required). The renormalization of the parton 
distributions is process independent, since the same collinear configurations occur in all 
processes. This is because they are related to hadronic emissions (or fragmentation into 
hadrons) which occur in a process independent manner and can be constructed from a uni- 
versal set of subdiagrams containing the configurations having collinear divergences with 
respect to the particular parton. The only difference then, from one process to the next, 
is the soft singularity structure. But the soft singularities cancel via the Bloch-Nordseick 
mechanism pSf (or KLN theorem [Q) such that the singular cross section for soft emissions 



is proportional to the Born term with (minus) the same overall factor as that coming from 
the loops. Hence, the scheme dependences of 0{e) which which multiply the 1/e singular 
terms, cancel exactly between the loops and the soft bremsstrahlung contributions. 

Thus, the only scheme dependence may come from the non-cancelling mass singularities, 
whose structure is process independent. Hence, all the argumentation here applies to all 
QCD processes at one-loop; extension to higher orders should be analogous. Also, the same 
conclusions concerning the regularization scheme independence in the MS^ scheme apply 
to all one-loop QCD processes. 
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V. DRED PARTON DISTRIBUTIONS AND 
FRAGMENTATION FUNCTIONS 



It is straightforward to show how to relate the DRED parton distributions and frag- 
mentation functions to those in DREG (MS). This is useful if one wishes to work strictly 
in DRED, but make use of the abundant sets of DREG parton distributions and fragmen- 
tation functions. First, we will consider the parton distributions, again dropping the A[t]'s 
for conciseness. 

We will make use of the fact that DRED is equivalent to the MS^ scheme of DREG. 
Noting that the bare parton distribution, /jy^(x), on the LHS of (^Oj) is factorization scheme 
independent, we obtain 

or 

= /S(^) = -Eg J'^flfAiy)m^/y) + ^(«^)- (79) 

From (^), we may immediately write for the fragmentation functions 

V^friz) = Vffliz) = Vf%z) - 5: g l'^V%iy)P!.{z/y) + 0{al). (80) 

j ^ y 

Comparing with ([76|) , ( [77|) we see explicitly that going from DRED to DREG MS sim- 
ply amounts to expressing the DRED (or MS^) parton distributions and fragmentation 
functions in terms of the DREG MS ones, and vice-versa from series inversion. 



VI. NUMERICAL RESULTS 



Here we present asymmetries and cross sections for the Drell-Yan process in p-p collisions 
at energies relevant to RHIC. In general, we use the two- loop MS expression for as(/i^), 
with four flavors and A = 0.2 GeV, except in the transversely polarized cross sections 
where we use the one-loop expression in order to be consistent with [§. Also, we take 
/i^ = Mj = M^. For the unpolarized cross sections, we use the DREG subprocess cross 
section convoluted with the unpolarized parton distributions of ||25| (Set S-MS). For the 
longitudinally polarized case, we use the MSg (or DRED) subprocess cross sections, since 
they are physically consistent (and regularization scheme independent), convoluted with the 
longitudinally polarized parton distributions of PU| (Set 1, SU{3) symmetric sea) which fit 
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well the recent DIS data except at very low x not covered for the kinematics considered 
here. 

For the transversely polarized subprocess cross sections, we again use the MS^ result. 
For the transversity distributions, we choose for the valence distributions (at Qq = 4 GeV^) 

AtF,„/p(x, Qg) = 2.lx-\l-xf-^ (81) 
^TFd^lp{x,Ql) = -.76x-«(l-x)3-4 (82) 

and for the sea-quark distributions 

^TFq/p{x,Ql) = -.l2x-\l-xf-\ q = u,d,s (83) 



(one half the value used in ^). These satisfy the upper bound proposed by Soffer 
As well, the valence distributions are consistent with bag model predictions. The small- 
and large-x behaviour is consistent with the longitudinal and unpolarized cases. There 
are no other definite constraints we may impose. Q^-dependent parameterizations for the 
longitudinal and transversity distributions are given in [^. We use one-loop evolutions, 
as two- loop polarized split functions do not yet exist. 
Define the asymmetries 



and 



da/dW^ 

Ma/dl 
da/dM^ 



Ada/dM^ , , 

= , , (85) 



noting that 

da 1 da 

dM^d(j)3 ^ 2^ZM2' ( ' 

Fig. 1(a) presents At for ^/S = 100 GeV, in the range 0.05 < y/r < 0.5. The largest 
value for —At is 8%. We notice that for the qq subprocess. At is reasonably perturbatively 
stable, unlike the cross sections which increase by 50 - 100 % under HOC [^]. Inclusion of 
the qg subprocess makes the asymmetry somewhat more negative since the qg subprocess 
contributes negatively as noticed in fl], |2|. 

Fig. 1(b) presents the corresponding polarized cross section, —A'^da/dM'^dip^. The 
sharp dropoff with increasing y/r is a combined result of the softness of Aq^Fg/p, the 1/M*^ 
behaviour of the cross section and the decreasing integration region with increasing ^/T. 

Fig. 2 presents the corresponding quantities for = 200 GeV and 0.05 < < 0.25 
(away from the Z-exchange region). Similar features hold, except the cross sections are 
somewhat smaller due to the suppression (which amounts to in A^^da / dMdip^). 

In Figures 3 and 4 we present the corresponding plots for Al and —Ada/dM"^. The 
largest value for —A^ is 16%. As expected from hehcity conservation, the qq subprocess 
exhibits great perturbative stability. Inclusion of the qg subprocess however, upsets this 
stability since it contributes with sign opposite to that of the qq subprocess (in both the 
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polarized and unpolarized cases) and is relatively large in the polarized case. Hence the 
net asymmetry becomes somewhat smaller in magnitude. This wouldn't be a problem in 
p-p collisions, where one is probing valence-valence distributions. For p-p collisions, the 
smallness of AFq makes the qg subprocess more significant. A smaller polarized gluon 
distribution and/or a larger polarized sea-quark distribution would reduce this effect. Still, 
for the larger ^/T, measurable asymmetries are obtained. 

For a discussion of the perturbative stability of the longitudinal asymmetry in p-p 
direct photon production, see |3^, ^ . There it was noted that the asymmetry is perturba- 
tively stable if the polarized gluon distribution is sufficiently large that the qg subprocess 
dominates. This contrasts with p-p Drell-Yan, where a large polarized gluon distribution 
destabilizes the asymmetry. 

VII. CONCLUSIONS 

We have presented complete NLO analytical results for the Drell-Yan process with 
unpolarized, longitudinally polarized and transversely polarized hadrons. These results are 
in a form valid for all n-dimensional schemes. It was shown how one can easily convert 
from results obtained in one scheme to those of another, regardless of the polarization, 
for one-loop QCD processes. This procedure simply amounts to expressing the parton 
distributions and fragmentation functions in one scheme in terms of those in the other. As 
well, the origin of the scheme dependences was elucidated. A mass factorization scheme, 
which we call MS^, was introduced. It was shown that in this factorization scheme, the 
final results are regularization scheme independent and coincide with those of DRED MS. 
A simple method for converting parton distributions and fragmentation functions from 
DREG to DRED was given. 

For p-p collisions at energies relevant to RHIC, asymmetries and cross sections for trans- 
versely and longitudinally polarized collisions were presented. For the transverse case, the 
asymmetries reached -8% and exhibited reasonable perturbative stablity. For the longitudi- 
nal case, the asymmetries reached -16% and the qq subprocess exhibited great perturbative 
stability. Inclusion of the qg subprocess somewhat lessened the longitudinal asymmetries 
however. Still, p-p collisions serve as the best probe of the polarized antiquark distributions 
in the proton, and they may be extracted with sufficient experimental statistics. 
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FIGURE CAPTIONS 

1. (a) Transverse asymmetry. At, in leading order and in next-to- leading order; (b) cor- 
responding next-to-leading order polarized cross section versus ^/T at \/S — 100 GeV. 

2. (a) Transverse asymmetry. At, in leading order and in next-to-leading order; (b) cor- 
responding next-to-leading order polarized cross section versus ^/r at = 200 GeV. 

3. (a) Longitudinal asymmetry, Al, in leading order and in next-to-leading order; 
(b) corresponding next-to-leading order polarized cross section versus at = 
100 GcV. 

4. (a) Longitudinal asymmetry, Al, in leading order and in next-to-leading order; 
(b) corresponding next-to-leading order polarized cross section versus a/t at y/S — 
200 GeV. 
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